Journal of Research of the National Bureau of Standards 



Vol. 54, No. 4, April 1955 



Research Paper 2583 



An Algorithm for Solving the Transportation Problem ' 

A. Gleyzal 2 

This article describes a new computational scheme for solving the transportation problem 
(described below) in which combinatorial ideas, rather than the theory of linear inequalities, 
play the major role. Whether the algorithm is superior to other methods previously 
employed or proposed must await computational experience, but preliminary inspection is 
encouraging. 



1. Introduction 

The transportation problem (stated for integers) 
is the following: Let C be an mXn matrix of integers 
en, and let a u . . . ,a m , b u . . ,,b n be positive 
integers such that a x -\- . . .+a m =bi + . . . + b n . 
We are required to find an mX% matrix X=(x ) 
such that 3 

(a) x tj is a nonnegative integer, 






Zu x ij — bj 



0=1, • . • ,n), 



and the linear form \ c {x)=^2 l c i jX i j is a minimum for 
. i-J ^'-**** i *** Mi 

all matrices satisfying (a) and (b) .1 1 

Definitions: A matrix X satisfying (a) and (b) is 
called a selection. If X and Y are selections, we say 
X(CHY(C) (read "X is preferable to Y with 
respect to C") if l c (x) ^l c (y)- If l c (z)<l c (2/), we say 
X(C)^Y(C)("X is strongly preferable to Y with 
respect to G Y "). It is clear that the relation ^ is 
transitive. If X is preferable to all selections with 
respect to C (i. e., X solves the problem), it is said to 
be optimal. 

If X is a selection, and i, j are indices such that 
Xij^l, then dj is called an X-selected element of C. 
If ri v . .,r m , s h . . .,s n are integers and D=(dij) = 
(Cij+rt+Sj), then C^D ("Cand D are equivalent")- 
It is easy to see that ~ is an equivalence relation. 
Further, one may prove that if C^D, then X(C)^ 
Y{C) if and only if X(D)<Y(D), and X{C) Y(Cj 



1 A code of this method in the case of the "personnel problem," using an 8X8 
matrix, has recently been composed for SWAC (National Bureau of Standards, 
Western Automatic Computer). For a variant of the method in this case, which 
is also applicable to matrices consisting of nonintegers, see Theodore S. Motzkin, 
The assignment problem (Proceedings of the American Mathematical Society 
Sixth Symposium for Applied Mathematics, 1954). 

2 Present address: Naval Ordnance Laboratory, White Oak, Silver Spring, 
Md. 

3 Other references to this problem and references to further literature are given 
by: George B. Dantzig, Maximimization of a linear function of variables sub- 
ject to linear inequalities (Simplex Method), chapter XXI, and Application of 
the Simplex Method to a transportation problem, chapter XXIII, both included 
in Activity analysis of production and allocation, edited by T. C. Koopmans 
(John Wiley & Sons, Inc., New York, N. Y., 1951); Merrill M. Flood, On the 
Hitchcock distribution problem, Pacific J. Math. 3, 369-386 (1953). 



if, and only if, X(D)^Y(D). Since throughout the 
computation subsequently described we work only 
with matrices equivalent to C, we shall henceforth 
write X(C)^Y(C) as X^Y (similarly, X<J). 

A closed circuit is a set of selected elements of the 
form c i]Jv c hJ2 ,c i2J2 ,c i2h ,. . .,c ith ,c UJV where i h . . .,i t 
are distinct indices and j h . . .,j t are distinct indices. 

2. Summary of the Computation 

The computation consists of two parts, which are 
performed alternately until an optimal selection is 
obtained. It will be seen that the algorithm itself 
yields the information that the optimal selection lias 
been reached. 

As a preliminary, we begin with some selection X. 
In part I, we find a matrix D^C and a selection 
Y^X such that the F-selected elements of D are 0. 
In part II, we discover if Y is optimal, and if not, 
obtain a selection Z-^Y. Putting Z in place of X, 
we begin again with part I and continue. As there 
are only a finite number of selections, it is clear from 
the fact that Z^(Y^X, that this process terminates 
in a finite number of steps, and we thus arrive at an 
optimal selection. 

3. Summary of Part I 

With the pair of matrices C, X, we associate a pair 
of nonnegative integers (c x ,x) defined as follows: 
c x is the number of nonzero X-selected elements of 
C, x is the number of nonzero elements of X. 

If c x =0 f we are finished with part I (set D=C, 
Y=X). If not, we construct, starting with a non- 
zero X-selected element of C a "tree I", which we 
discover to be either in case I A or case IB. If in 
case I A, we find a matrix C such that C'~C and the 
number of nonzero X-selected elements of C f is less 
than c x . If in case IB, we find a selection X' 
preferable to X such that the set of nonzero elements 
of X f is a proper subset of the set of nonzero elements 
of X. Thus, after completing I A or IB, we have a 
matrix C ' ^C and a selection X'^X such that 



c'z'^c x , x f ^x, 



a) 
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and strict inequality holds in at least one of the two. 
If <v>0, we begin again, obtaining C" , X" ' , and 
so on. By virtue of (1), this cannot continue 
indefinitely. But we cannot have & (n) = 0, for this is 
inconsistent with (b); therefore, for some n, c xin) =0. 
Hence, the desired result of part I is achieved. 

4. Summary of Part II 

In part II, we begin with a pair of matrices D, Y 
such that d y =0. Let e D be the minimum of all 
elements of D that are not F-selected. Let f D be 
the number of elements of D equal to e D . 

If e D ztO, then, clearly, Y is optimal. Assume the 
contrary, therefore. Starting with an element 
d tiJl of D such that d tl j =e D) we construct a "tree II." 
We examine this tree and discover whether we are 
in case HA or case IIB. If in case HA, we construct 
a matrix D'~D such that the F-selected elements 
of D' are 0, e D ^e D (and the behavior of f D is un- 
specified), or 



e D > = e D , Jd'KJd, 



(2) 



and we put D' in place of D and begin again w4th 
part II. It is obvious from (2) that we cannot 
continue to be in case HA indefinitely, for eventually 
we would obtain e D (n)^0, which implies that Y is 
optimal. Therefore, if Y is not optimal, we must 
at some stage enter case IIB. If in case IIB, we 
find a selection Z-^Y, and go back to part I, with 
any matrix equivalent to C (C itself, if convenient) 
replacing C, and Z replacing X. 

5. Detail of Part I 

Assume c ioJQ is an -XT-selected element of C that is 
not zero. We define, inductively, the following 
"family tree I": 

1. The "founder" is c iQJ0 ("first generation"). 

2. Assume now that we have defined the kth 
generation, for Jt=l, . . .,«. We now define the 
(s+l)th generation; i. e., the union, over all mem- 
bers of the sth. generation of the "children" of each 
member of the sth. generation. The "children" of 
such a member are: if s is odd, all other X-selected 
elements in the same row of C as the given member; 
if s is even, all other X-selected elements in the same 
column of C as the given member. If the (s+l)th 
generation is vacuous, the tree terminates with the 
5th generation, and we say we are in case IA; if the 
(s+l)th generation contains at least one member of 
a previous generation (possibly the sth) the tree 
terminates with the (s + l)th generation, and we say 
we are in case IB. If the tree does not terminate, 
we continue. The tree must terminate, of course, 
after a finite number of generations. 

Case IA: In this case, we shall show how to con- 
struct numbers r ly . . ,,r m , $i, , . .,s n such that 
C / ={c f ii )={c i j J rr i -\-s j )^C y and c x <c % (indeed, all 
the elements of C corresponding to the elements of 
the tree will be zero). 



Associated with each element of the tree, define 
inductively the following numbers N(t): 

If p is the parent of t, define 

N(t) = -t-N(p). 

We note that because we are in I A, rather than IB, 
every element has a unique parent. Now, let a 
be any row index of some element c a p of an odd 
generation: for each such a, define r a =N(c a p). 
Observe that because we are in case IA, this defini- 
tion is unambiguous. For suppose c a ^ and c a ^ 2 were 
two elements of row a, each occurring in an odd 
generation, then it follows from the definition of the 
tree that since c aj g occurs in an odd generation, c a p 2 
is its child. Similarly, c a pi is the child of c a p . This 
would mean that we were in case IB. 

For every index i=l, . . .,m, which is not a row 
index of an element of an odd generation, define 
r 4 =0. Similarly, let /3 be any column index appear- 
ing in some element c a p of an even generation: for 
each such /3, define sp=N(c a p). For each index 
j=l, . . .,n, which is not a column index of an 
element of an even generation, define fy=0. 

It is now easily seen that C has the desired 
property. 

Case IB: In this case, we shall show how to find 
a selection X'^KX such that x'<Cx. We first show 
how to find among the members of the tree a closed 
circuit. Let k be defined by the statement that the 
(& + l)th generation is the first one that contains an 
element e of the rth generation, where r<^k + l, and 
let e k be a member of the kth. generation whose child 
is e T . (Although the knowledge is not necessary for 
the subsequent argument, it is convenient to note 
that & + l>3.) In other words, it is at the (& + l)th 
generation that we discover that we are in case IB. 
Hence each member of generations k, k—1, ... ,2 
has a unique parent in the preceding generation, for 
if a member of the tree has two parents in the pre- 
ceding generation, then all three members are in the 
same row or column, one of the parents would be a 
child of the other parent, and we would have arrived 
at case IB earlier. 

Now let e k _i be the unique parent in the (k — l)th 
generation of e k , e k - 2 be the unique parent in the 
(k—2)th generation of e k _ h etc. There are two 
possibilities: either (a) e r is an ancestor of e k , or (b) 
e r is not an ancestor of e k . If (a) holds, then clearly 
e k} e k -i, e k _ 2 , • • • A+b e r is a closed circuit. If (b) 
holds, let s be defined by the statement that the sth 
generation is the last generation that contains a 
common ancestor e 8 of both e k and e r . Let e' r _ 1 be 
the parent of e r in the (r— l)th generation, e'_ 2 the 
parent of e,_ x , and so on until we reach e 8 . From 
this "line of descent" e r) e r _ u 6 r '_ 2 , . . . , e' s+1 , e s and 
the previous e k , e k ^ 1} e k _ 2 , • • • , e«-H> e 8 , noting that 
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e' a+ i, e 8 , e s+1 are in the same row or column, it is easy 
to see that 



e ki €k-l> € k-2y • • • j ^s+l> ?s 



@r—2) &r-l) &r 



is a closed circuit. For the definition of k implies 
that the same row or column index occurs only in 
consecutive elements of this cycle, and if we con- 
sider these elements as arranged in a cyclic order, 
the construction of the tree and definition of k 
imply that the same row or column index occurs 
only in consecutive elements. 

Let us denote this closed circuit by 



c h*V Ci ih 



$i t h> 



as explained in the introductory definitions. 
Let 



^2=c ilh +e hh + .,.+€, 



/•V 



Either 2i^S 2 or 2!>S 2 . Assume the former, and 
let % 2 =mm(x hhJ . . . ,a? Vl ). 
Con side i- the new selection X' 

For all other pairs of indices i, j, #»•/=#,•/. 

It is easy to verify that X' is a selection (i. e., 
satisfies (a) and (b), that X'^sX, and #'<#. 

If 2!>2 2 , let x 1 =mm(x iiJ - i , . . . ,aP^ t ). Consider 
the new selection X' . 

*i'i = "^*i'i ^ lf ' ' ' > ^ { t 1 ' t = ^ { t ] ' t ^ 1 

^t 1 7 2 ==^t 1 7 2 + ^li • • • , #• !/! = #»,>*! + #!• 

For all other pairs of indices i, j, £•,•=:*;,•. Then it is 
easy to verify that X' is a selection, X'-^X, a;'<jr. 
This completes the discussion of part I. 

6. Detail of Part II 

Assume d ilfl is an element of D that is not F- 
selected, and such that 0>^ i7 - =e D = minimum of all 
elements of D not selected by F. We construct in- 
ductively, the following "family tree II". 

1. The founder is d i{ji . 

2. Assume now that we have defined the kth 
generation, for k=l, . . . ,s. We now define the 
(s + l)th generation, the union, over all members of 
the sth generation of the children of each member of 
the sth generation. These are the children: If s 
is odd, s^>l and there exists a F-selected element of 
D whose column index is j u then this element is its 
child, the tree terminates, and we say we are in case 
IIB. If no such element exists or if 8=1, then the 
children consist of all F-selected elements in the 



same row as the parent member, omitting, however, 
those whose column index is the same as that of a 
member of generations 1, . . . ,s. If s is even, the 
children of a member are all aonpositive elements of 
D in the same column as the member, omitting, how- 
ever, those whose row index is the same as that of a 
member of generations of 1, . . . ,s. 

Eventually, the tree must terminate, and if it does 
not terminate in case IIB, we say we are in case II A. 

Case IIA: In this case, we shall show how to 
construct numbers r u . . . ,r w , s 1} . . . , s n such 
that D , = (di i ) = (d ij + r i +s i )^D, and D' has the 
properties (2). 

Let a be any row index appearing in some number 
d a p of an odd generation. For all such a, define 
r a = J rl. For any index i=l, . . . , m not in- 
cluded in [a], define ^=0. Similarly, let P be any 
column index appearing on some member d a $ of an 
even generation. For all such 13, define s j8 = — 1. 
For any index j—l, . . . , n not included in {/3}, 
define s ; = 0. 

We now show that I)' has the properties (2). 
First, the F-selected elements of IY are 0: For if 
(Iy is such an element, and if i resp. j appears as a 
row resp. column index of a member of an odd 
resp. even generation, then j resp. i appears as a 
column resp. row index of a, member of an even 
resp. odd generation; hence, ^+^=0. Second, if 
dij<^0 and s ; = — 1, then by our construction r<=l; 
hence the minimum of the elements of D f is not 
less than the minimum of the elements of D. In ad- 
dition, because we are in case IIA, d' { ,• =d i{i -\-l^>e D . 
These two statements imply the remainder of (2). 

Case IIB: Because we are in case IIB, it is clear 
how to find a closed circuit in the tree containing 
d t * . Let us denote this circuit by 

Then d^d^ , . . •Ai t j l are F-selected elements of 
D. Consider the new selection Z: 

For all other pairs of indices i,j, define 2ij=yij. 

Then it is easy to verify that Z is a selection. 
Further, because all elements of the circuit are non- 
positive, d iih = • • • =d itJi = 0, and d Vi <0, it 
follows that Z^X. 

7. An Example 

Consider the transportation problem with 



(7= 



"2595^ 
8358 
7314 
5972J 



and a,= (3,2,3,3), ^=(3,5,2,1). For purposes of 
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illustration, we assume an initial selection matrix, 



X= 



nioo^ 

0020 
0300 

Uioo 



It is convenient to indicate the selection by dots 
placed above the elements of the matrix C, as follows, 



f 2 5 9 5" 

8 3 5 8 



7 3 14 
L5 9 7 2J 



The steps in the solution and the operations associ- 
ated with cases IA, IB, IIA, and IIB will be clear 
from the following: 4 



The last matrix is a solution. No closed circuit of 
zeros exists. Hence the solution is unique. For pur- 
poses of illustration, we have followed the method 
explicitly. If one permits some flexibility in choosing 
the sequence of operations, much faster convergence 
may be obtained. 



The author expresses his gratitude to Alan J. Hoff- 
man to whom he is indebted for transforming the 
method and proof into the algebraic form given here. 
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(?) 5 9 5 




8 3 5 8 




5 9 7 2 







3 7 3 




3-203 


-3 


7 © 1 4 




5 9 7 2 






-1 
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-i 
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-I 
9 


4 

/ 


r2~ 


-^6 


10 


■( 


f 


-2 


8 


©_ 


\ 


-2 


6 





* Numbers on borders at matrices are quantities added to the indicated row 



'.Numbers on borders at matrices arc Quantities added to tne indicated row ttt rjuTTvinmnAT At^tmI 9Q IQ^A 

column. Circles drawn about an element indicate the first element of the tree. WASHINGTON, April Zo, iy04t. 
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